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Avian coronavirus infectious bronchitis virus (IBV) is a respiratory pathogen of chickens, causing severe 
economic losses in poultry industry worldwide. Live attenuated viruses are widely used in both the broi¬ 
ler and layer industry because of their efficacy and ability to be mass applied. Recently, we established a 
novel reverse genetics system based on targeted RNA recombination to manipulate the genome of IBV 
strain H52. Here we explore the possibilities to attenuate IBV in a rational way in order to generate safe 
and effective vaccines against virulent IBV (van Beurden et al„ 2017). To this end, we deleted the 
nonessential group-specific accessory genes 3 and/or 5 in the IBV genome by targeted RNA recombination 
and selected the recombinant viruses in embryonated eggs. The resulting recombinant (r) rIBV-A3ab, 
rIBV-A5ab, and rIBV-A3ab5ab could be rescued and grew to the same virus titer as recombinant and wild 
type IBV strain H52. Thus, genes 3ab and 5ab are not essential for replication in ovo. When administered 
to one-day-old chickens, rIBV-A3ab, rIBV-A5ab, and rIBV-A3ab5ab showed reduced ciliostasis as com¬ 
pared to rIBV H52 and wild type H52, indicating that the accessory genes contribute to the pathogenicity 
of IBV. After homologous challenge with the virulent IBV strain M41, all vaccinated chickens were pro¬ 
tected against disease based on reduced loss of ciliary movement in the trachea compared to the non- 
vaccinated but challenged controls. Taken together, deletion of accessory genes 3ab and/or 5ab in IBV 
resulted in mutant viruses with an attenuated phenotype and the ability to induce protection in chickens. 
Hence, targeted RNA recombination based on virulent IBV provides opportunities for the development of 
a next generation of rationally designed live attenuated IBV vaccines. 

© 2018 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license 

(http://creativecommons.Org/licenses/by-nc-nd/4.0/). 


1. Introduction 

Infectious bronchitis virus (IBV) is an avian gammacoronavirus 
that belongs to the family Coronaviridae of the order Nidovirales 
[1,2], It was first discovered in the United States in the 1930s [3] 
as the causative agent of a highly contagious respiratory disease 
in chickens, known as infectious bronchitis. Although IBV princi¬ 
pally infects the upper respiratory tract, some IBV strains affect 
the renal tubuli, oviduct and parts of the gastrointestinal tract. 
Infection with IBV may lead to reduced growth and egg production, 
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and is as such regarded as one of the economically most relevant 
viral pathogens in the poultry industry worldwide. 

Current strategies to prevent IBV in poultry include vaccination 
with live attenuated IBV vaccines as well as with inactivated vac¬ 
cines. The widely used live attenuated IBV vaccine, HI20, was 
developed in the 1960s in The Netherlands by serial passaging of 
a Massachusetts-like IBV strain in embryonated eggs [4], After 
120 passages, the resulting virus H120 had become strongly atten¬ 
uated as a result of its embryo adaptation and did not cause signif¬ 
icant disease in young chicks, while it induced an immune 
response protective against challenge with wild type Mas¬ 
sachusetts IBV [4], The H120 vaccine has been used successfully 
for decades. With the worldwide occurrence of IBV in both com¬ 
mercial and backyard chicken in a wide variety of geno-, sero- 
and protectotypes [5], protection against IBV by vaccination has 
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become far more complicated nowadays. Yet, the development of 
new live attenuated vaccines is still done by serial passaging of 
virus field isolates in embryonated eggs, which is a laborious and 
time-consuming process with unpredictable outcome with regards 
to vaccine safety and attenuation. 

In order to investigate functions of 1BV proteins in a directed 
way, several research groups independently developed systems 
to manipulate the IBV genome 6-10]. These reverse genetics sys¬ 
tems (RGS) are based either on the cell-culture adapted and highly 
attenuated IBV strain Beaudette or on the attenuated IBV vaccine 
strain H120. Since IBV Beaudette and H120 are highly attenuated 
[4,11 ], the ability to study virulence factors of IBV in vivo is limited. 
This would require the introduction, or substitution, of factors con¬ 
tributing to the infection in vivo, including (parts of) the replicase 
genes and for example spikes from other IBV serotypes [12-14], 
We have recently established a reverse genetics system based on 
the more virulent IBV strain H52 [4,15]. Using targeted RNA recom¬ 
bination and performing the rescue and selection of candidate 
recombinants in embryonated chicken eggs, we solved the bottle¬ 
neck of the inability to propagate virulent IBV strains to grow in 
cell culture. This provided a novel way to generate recombinant 
IBV for in vivo studies, including vaccine development. 

The IBV genome encodes the nonstructural proteins involved in 
replication of the viral genome (ORFlab) and the structural pro¬ 
teins spike S, membrane M, envelope E, and nucleocapsid N. In 
addition, genes 3 and 5, located between S and M and M and N, 
respectively, code for the proteins 3a, 3b, and 3c (E), and 5a and 
5b [2], These gammacoronavirus-specific proteins 3ab and 5ab 
are non-structural and non-essential for virus replication 
[9,16,17], Several lines of evidence indicate, however, that they 
might have a role in viral pathogenesis in chickens. In particular, 
these genes are conserved across all IBV field strains [18] and 
recent studies indicate that the 3a and 3b proteins induced a 
delayed activation of the type I interferon (IFN) response in vitro, 
with protein 3a additionally being involved in resistance of IBV 
to the cellular antiviral state induced by IFN [19,20 . Accessory pro¬ 
tein 5b was found to contribute to host cell shut-off, including 
amongst others the inhibition of translation of type 1 IFN 21 ]. 
Based on these observations, and on the reported role of the acces¬ 
sory genes of other coronaviruses in pathogenicity [22-25], we 
reasoned that the deletion of the accessory genes 3 and 5 might 
attenuate the more virulent phenotype of IBV H52 in vivo. Here 
our RGS [15] was used to generate mutant recombinant IBVs for 
vaccine development based on IBV H52 BI. We show that recombi¬ 
nant viruses lacking the 3ab and/or the 5ab gene cluster are viable 
and replicate like HBV-H52, while showing reduced pathogenicity 
in chickens. Their ability to protect vaccinated chickens against a 
homologous challenge demonstrates the feasibility of this 
approach in creating next generation vaccines against infectious 
bronchitis. 


2. Materials and methods 

2.1. Cells, eggs & viruses 

Murine LR7 cells [26] were cultured in Dulbecco’s Modified 
Eagle Medium (DMEM; BioWhittaker), supplemented with 4 mM 
1-glutamine (Lonza), 10% Fetal Calf Serum (FCS; BioWhittaker) 
and 0.05 mg/ml gentamicin (Gibco Invitrogen), at 37.0 °C and 5% 
C0 2 . 

Fertilized specific pathogen free (SPF) white leghorn eggs (Ani¬ 
mal Health Service, Deventer, The Netherlands) were incubated at 
37.5 °C and 45-65% relative humidity. Embryonated chicken eggs 
(ECE) were inoculated into the allantoic cavity at day ten of incu¬ 
bation, unless stated otherwise, and candled twice daily. Eggs were 


transferred to 4 °C for 16-24 h prior to collection of the allantoic 
fluid (AF) and the chorio-allantoic membrane (CAM). Virus titra¬ 
tion in ovo was based on 50% embryonic infectious dose (E1D 50 ) 
per ml, as determined at day 7 post inoculation (pi) according to 
Reed and Muench 27], For preparing virus stocks, AF of four to 
ten ECEs inoculated with 100 EID 50 were pooled after incubation 
for 24 h. 

IBV strain H52 BI (Boehringer Ingelheim, Bl, Ingelheim, Ger¬ 
many), recombinant IBV wild-type derived from H52 (rIBV-wt) 
and murinized (m)IBV (strain #1B3-I1A) were propagated as 
described before [15], 

2.2. Immunohistochemistry (IHC) 

Immunohistochemistry on CAMs collected from ECEs was per¬ 
formed as described previously 15] using monoclonal antibody 
(MAb) Ch/IBV 26.1 directed against the IBV S2 protein (Prionics, 
Lelystad, The Netherlands) [28,29], 

2.3. RNA isolation, reverse transcription & PCR 

RNA was isolated using the QIAamp viral RNA Mini Kit (Qiagen, 
Hilden, Germany) according to manufacturer’s protocol. Reverse 
transcription was performed using the Transcriptor First Strand 
cDNA Synthesis Kit (Roche, Basel, Switzerland) according to manu¬ 
facturer’s protocol, with random hexamers. PCR was performed 
with Phusion Hot Start II High-Fidelity DNA Polymerase (Thermo 
Fisher Scientific) for sequencing and cloning purposes. One-step 
RT-qPCR was used to semi-quantitatively assess virus load in AF, 
using primers IBV.RdRp.F41 and IBV.RdRp.R41 as previously 
described [15]. 

2.4. Construction of p-lBV-A3ab, p-IBV-A5ab, and p-lBV-A3ab5ab 

The design of donor plasmid p-IBV has been described previ¬ 
ously [15] (Fig. 1A). The constructs p-IBV-A3ab, p-IBV-A5ab, and 
p-IBV-A3ab5ab, in which the accessory genes 3ab, 5ab, and 3ab5ab 
are deleted, are shown in Fig. IB and C. Design of the delta3ab frag¬ 
ment (A3ab) was such that the 1 nt overlap between the stop 
codon of the spike gene and the start codon of the 3a gene was 
replaced by an overlap between the stop codon of the spike gene 
and the start codon of the E gene (Fig. 1C-5). For the delta5ab frag¬ 
ment (A5ab) the start codon of the 5a gene is now used as start 
codon of the nucleocapsid gene (Fig. 1C-6). DNA fragments span¬ 
ning the semi-unique restriction enzyme sites (RES) surrounding 
the accessory genes with the deletions designed as described 
(A3ab and A5ab) were cloned into pUC57-simple by Genscript 
(Piscataway, NJ, USA). Delta3ab was ligated into p-IBV-5-lb-S-SIR 
[15 after Nhel-Pmll double digestion to remove the 3ab gene, fol¬ 
lowed by the 3T region to form p-IBV-5-lb-S-SIRA3ab-3T, now 
called p-IBV-A3ab. Likewise, A5ab was ligated into p-IBV-3T [15] 
after Afel-Nhel double digestion to remove the 5ab gene. Fragment 
3TA5ab was subsequently ligated after p-IBV-5-lb-S-SIR or p-IBV- 
5-lb-S-SIRA3ab to form p-IBV-5-lb-S-SIR-3TA5ab or p-lBV-5-lb- 
S-SIRA3ab-3TA5ab, respectively, now called p-IBV-A5ab and p- 
IBV-A3ab5ab. Composition of each of the plasmids was confirmed 
by PCR, restriction enzyme digestion and sequencing of the inserts 
(Macrogen, Amsterdam, The Netherlands) (see Table 1). 

2.5. Targeted RNA recombination and rescue of recombinant IBVs 

rIBVs were generated by introducing the IBV spike ectodomain 
into the mIBV genome by targeted RNA recombination between p- 
IBV donor RNAs and recipient virus mIBV, as described previously 
[15]. In short, capped run-off donor transcripts were synthesized 
from p-IBV-A3ab, p-IBV-A5ab and p-IBV-A3ab5ab using the 
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O Spike endodomain - 3a junction 

IBV CTGAACAATACAGACCTAAAAAGTCTGTTTAATGATCC 

pIBV CTGAACAATACAGACCTAAAAAGTCTGTTTAATGATCC 

Spike TEQYRPKKSV . 


3a 


M I 


© 3b - Envelope junction 
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p-IBV 


CAGGTG liATGAATTTATTGAATAAGTCGCT 
CAGGTGATGATGAATTTATTGAATAAGTCGCT 


3b 

Envelope 
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M M N L L N 


© Intergenic region - 5a junction 

IBV CTTAACAAATACAGACG ATGAAATGG 

p-IBV CTTAACAAATACAGACGAT AAATGG 

5a M K W 


© 5b - Nucleocapsid junction 

IBV TTTGTCA TG GCGAGC [...] ATCAAACTAGGAGGGC 

p-IBV TTTGTCATGGCGAGC [...] ATCAAACTAGGAGGGC 

5b LSWRA SN. 

Nucleocapsid MAS IKLGG 

e A3ab junction 

p-IBV-A3ab 
Spike 
(3b) 

Envelope 

Q A5ab junction 

p-IBV-A5ab CTTAACAAATACAGACG ATGGCGAGC [...] TCAAACTAGGAGGGC 

Nucleocapsid MAS IKLGG 

(5b) W R A S N . 


CTGAACAATACAGACCTAAAAAGTCTGTTTA A^ i ATGAATTTATTGAATAAGTCGCT 
TEQYRPKKSV . 

N D E F I E 
MMNLLNKSL 


Fig. 1. Schematic overview of targeted RNA recombination principle and donor plasmids. (A) Schematic overview of step 2 in the targeted RNA recombination method for 
generating recombinant (r)IBV wild type (wt) [15]. IBV sequences are represented in blue, MHV sequences in red. pIBV-derived synthetic donor RNA is indicated by a black 
composite line above which the parts derived from specific sub-plasmids are indicated. PCR amplicons used to confirm the recombination (set C) and the status of gene 3ab 
(set E) and gene 5ab (set I) are depicted as black bars drawn to scale above the rIBV-wt genome, with encircled letters referring to the primer sets in Table 2. (B) Schematic 
layout of the donor plasmids p-IBV-A3ab, -A5ab and -A3ab5ab used in targeted RNA recombination to generate rIBV-A3ab, -A5ab and -A3ab5ab, respectively. (C) Nucleotide 
sequences of the gene and plasmid junctions are marked with numbers corresponding to the numbers in black circles in the schematic donor plasmid layout in (A) and (B). 
Transcription regulatory sequences are in bold. Nucleotide sequences are indicated for wild type IBV and donor plasmids p-IBV, p-IBV-A3ab and p-IBV-A5ab. Start codons are 
highlighted in green, stop codons in red, and start-stop overlaps in yellow. ORFs are underlined, overlapping ORFs are underlined with a bold line, and ORF translations are 
indicated as amino acids below the nucleotide sequences if applicable. (For interpretation of the references to colour in this figure legend, the reader is referred to the web 
version of this article.) 
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Table 1 

Plasmids used for generation of the donor plasmids p-IBV-A3ab, -A5ab and -A3ab5ab. 


Plasmid 

Genes 

Coordinates 

Length (nt) 

3'-end RES 

Surrounding RES 

Inserted in plasmid 

p-IBV-5 

5'-UTR 

1...497 

497 

BstBI 

n.a. 

n.a. 

p-IBV-lb 

lb, S 

19,610...20,379 

770 

Xhol 

BstBI 

p-IBV-5 

p-IBV-S 

s 

20,379.. .23,590 

3211 

Sty\ 

Xhol 

p-IBV-5-lb 

p-lBV-SIR 

S, 3ab, E, M 

23,591...25,318 

1728 

EcoRl 

Styl 

p-IBV-5-lb-S 

p-lBV-3 T 

5ab, N, 3'-UTR 

25,319...27,730 

2422 

MssI, Pac\ 

EcoRl 

p-IBV-5-lb-S-SIR 

p-A3ab 

n.a. 

23,604.. .24,304 

349 

n.a. 

Nhe I, Pmll 

p-IBV-5-lb-S-SIR 

p-A5ab 

n.a. 

25,467.. .25,974 

115 

n.a. 

Afe I, Nhe I 

p-IBV-3 T 


mMessage mMachine T7 kit (Ambion by Thermo Fisher Scientific) 
after Mssl-linearization. In vitro transcribed RNA was transfected 
by electroporation into LR7 cells previously infected with mIBV fol¬ 
lowed by intra-allantoic injection in ten-day-old embryonated SPF 
chicken eggs using five eggs per dilution (10 _1 -10 -5 ). Selection of 
candidates was done based on IHC and RT-qPCR. Suitable candi¬ 
date recombinants were subjected to two additional rounds of 
end-point dilution in eight-day-old ECE, and virus stocks were 
grown and titrated. The nucleotide sequence of the genomes of 
the rIBV delta variants region starting from nt 18,612 onwards 
was confirmed by RT-PCR and Sanger sequencing, using the primer 
sets specified in Table 2. 

2.6. In ovo growth kinetics of rlBV-A3ab, rIBV-A5ab, and rlBV- 
A3ab5ab 

The growth of i!BV-A3ab, iTBV-A5ab and riBV-A3ab5ab was 
compared to that of IBV H52 BI and rIBV-wt by performing RT- 
qPCR on viral RNA extracted from AF of inoculated eight-day-old 
ECEs. A 10-fold dilution series of IBV H52 BI RNA was used as ref¬ 
erence for quantification of EID 50 /ml equivalents, as described pre¬ 
viously [15], 

2.7. In vivo vaccination-challenge experiment 

A vaccination-challenge experiment was performed in SPF 
layer-type chickens, housed in isolators under controlled condi¬ 
tions, including HEPA-filtered supply and exhaust air. One-day- 
old chickens were inoculated via eye-drop with 10 3 EID 50 in 0.1 
ml of either wild-type IBV H52 BI (n= 15) or recombinant IBVs 
(four groups of n = 15 each). As controls, a non-vaccinated but 
challenged group (n = 10) and a negative, non-vaccinated, non- 
challenged, control group (n = 10) were included. The results from 
IBV H52 BI and rIBV-wt have been reported before [15]. Seven days 
post vaccination, five animals were removed from each of the five 
vaccinated groups and the negative control group. After euthanasia 
the tracheas were evaluated for ciliary activity as described previ¬ 
ously [15], 

After 21 days, the remaining animals of the five vaccinated 
groups (n = 10 each) and the challenge control group (n = 10) were 


challenged via eye-drop with 10 3 EID 50 in 0.1 ml of the virulent IBV 
strain M41 (Animal Health Service, Deventer, The Netherlands). 
Seven days post challenge, all remaining animals were euthanized, 
and their tracheal ciliary activity was evaluated. A one-way ANOVA 
with Tukey’s honestly significant difference (HSD) post hoc test 
was performed to analyze whether the ciliostasis scores between 
the groups were statistically different. 

Analysis of the recombinant viruses was done based on ciliosta¬ 
sis assays as recommended for live attenuated IBV vaccines by the 
European Pharmacopoeia. During the course of the experiment, 
birds were inspected on general and specific respiratory clinical 
signs, including gasping, coughing, tracheal rales, and nasal 
discharge. 

3. Results 

3.1. Generation of rIBV-A3ab, rIBV-A5ab, and rIBV-A3ab5ab 

LR7 cells infected with mIBV and transfected with in vitro tran¬ 
scribed donor RNA from p-IBV-A3ab, p-IBV-A5ab, and p-IBV- 
A3ab5ab were inoculated into the allantoic cavity of ten-day-old 
ECEs. At seven days p.i., with the exception of one embryo in the 
IBV-A3ab group, no embryonic death was observed, but several 
embryos in groups inoculated with the lowest dilutions showed 
stunting and curling typical for IBV. The presence of replicating 
IBV was investigated by IHC on CAM tissues using an antibody 
against the IBV S2 protein (Fig. 2). CAMs from eggs that received 
non-electroporated control cells showed no viral antigen produc¬ 
tion (Fig. 2, bottom row) while inoculation with infected and trans¬ 
fected cells in all cases resulted in the production of viral antigens 
in the CAM as detected by immunohistochemistry (Fig. 2, upper 
row). The presence of IBV RNA in the AF of these eggs was con¬ 
firmed by RT-qPCR (data not shown). 

3.2. Genetic characterization of r!BV-A3ab, rIBV-A5ab, and rlBV- 
A3ab5ab 

The genetic identity of rlBV-A3ab, rIBV-A5ab, and rlBV- 
A3ab5ab was assessed by RT-PCR using the primers specified in 
Fig. 1 and Table 2. Primer set [C] targeted the lb-spike gene junc- 


Table 2 

Primer sets used for the characterization of the 3' 9 kb of the viral genome of rIBVs. 


Primer set 

Target 

Primer 

Sequence (5' -> 3') 

Amplicon (bp) a 

C 

IBV lb - IBV S 

IBV.F73 

IBV.R73 

TCAGCATGGACGTGTGGTTA 

CCCCATGTAAATGCCAACCA 

992 

E 

IBV S - IBV M 

IBV.F15 

IBV.R15 

TGCTGCTTCCTTTAATAAG 

CTGCGACAAGACCTCCTG 

1994/1646 

I 

IBV IR - IBV N 

IBV.F28 

IBV.R16 

TGTTGTAGGTTGTGGTCCCA 

CTGAGGTCAATGCCTTATC 

1633/1248 


a For primer sets [E] and [I], the amplicon lengths for rIBV with/without gene 3 and/or 5 are given. 
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rIBV-wt rlBV-A3ab rlBV-A5ab rlBV-A5ab 



Fig. 2. Immunohistochemistry of CAMs after rescue of rIBV-A3ab, rIBV-A5ab and rlBV-A3ab5ab. Embryonated chicken eggs were inoculated with mIBV-infected LR7 cells 
that had been transfected with transcripts from donor plasmids p-IBV-A3ab, -A5ab or -A3ab5ab by electroporation. mIBV-infected non-transfected LR7 cells served as 
controls. Formalin-fixed and paraffin-embedded CAMs were stained using a monoclonal antibody against 1BV-S2. 


tion, using a forward primer located upstream of the BsfBI RES in 
gene lb, with the aim of confirming the correct recombination 
between mlBV and the donor RNA; primer set [E] spanned gene 
3 (resulting in a 348 nt shorter amplicon in case of 3ab deletion); 
primer set [1] spanned gene 5 (resulting in a 385 nt shorter ampli¬ 
con in case of 5ab deletion). The RT-PCR on viral RNA extracted 
from AF from ECEs inoculated with the infected and transfected 
LR7 cells showed for rIBV-A5ab and rlBV-A3ab5ab, but not for 
rlBV-A3ab, amplicons of different sizes across the mutated regions 
(data not shown). AF from eggs with a single RT-PCR product of the 
expected sizes (Fig. 3, passage (p) 0) were selected for further use. 


3.3. rIBV-A3ab, rIBV-A5ab, and rlBV-A3ab5ab are genetically stable 

Endpoint dilutions were performed starting with the initial AF 
stock that was positive in 1HC and RT-qPCR. The virus stocks were 
passaged two additional times on ECEs to generate a virus seed and 
working stock. All four passages of rlBV-wt, rlBV-A3ab, i!BV-A5ab, 
and i!BV-A3ab5ab were analyzed by RT-PCR using the primer sets 
[C], [E], and [I] (Fig. 3). No additional PCR fragments or fragments of 
other sizes were observed during passaging, confirming that the 
deletions were stably maintained in embryonated eggs across 
these passages. 



Fig. 3. Genetic characterization of rIBVs. PCR was performed on plasmid DNA (p) and on cDNA templates of viral RNA extracted from AF of ECEs inoculated with mIBV- 
infected and donor plasmid transcript-transfected LR7 cells (0), and subsequent passages (1-4) of rIBV-wt, rIBV-A3ab, rlBV-A5ab and rIBV-A3ab5ab. Primer set letters 
correspond to the respective locations in Fig. 1, and include primer sets that span the genomic area of recombination (C), accessory genes 3 (E) and 5 (I). Expected amplicon 
sizes are indicated in the right column, indicating for primer sets E and I the length in the presence and absence of genes 3ab or 5ab, respectively. M: molecular weight 
marker. 
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Sequence analysis of the 3' 9 kb of rIBV-A3ab, rIBV-A5ab, and 
rlBV-A3ab5ab (starting 1 kb before the BstBI RES in gene lb) con¬ 
firmed the intended genomic sequence of the viruses in the 
working stock (Supplementary file SI and Fig. 1C). One non- 
synonymous (A/C) substitution was observed in the spike of 
rIBV-A5ab P4 at position 23,303 (Ala -> Glu). At position 25,461 
in the intergenic region of rIBV-A5ab a C to A substitution 
was apparent. In rlBV-A3ab5ab, a non-synonymous A to G substi¬ 
tution was observed in the spike gene at position 22,076 
(Asn -> Ser). 


3.4. Recombinant IBVs have comparable growth kinetics in 
embryonated eggs 

The in ovo growth kinetics of the rIBVs were assessed by inocu¬ 
lating ECEs with 10 2 EID 50 of working stock per egg, and determin¬ 
ing the relative viral load in the AF of five eggs per virus strain at 6- 
12 h intervals by RT-qPCR. At the earlier time points 6 and 12 hpi 
the viral loads observed for rIBV-A5ab and rIBV-A3ab5ab were 
somewhat lower than for IBV H52 BI, rIBV-wt and rIBV-A3ab 
(Fig. 4A). At 48 hpi, the time when most embryos had died 
(Fig. 4B), viral loads were comparable for all viruses tested, indicat¬ 
ing that the viruses propagated to similar titers. 




0 6 12 24 36 48 

Time p.i. (h) 


Fig. 4. In ovo characteristics of rIBVs. (A) Quantitative RT-qPCR analysis of RNA 
extracted from AF of ECEs collected at 12 h intervals after inoculation with IBV H52 
BI, rIBV-wt, rIBV-A3ab, rIBV-A5ab, and rIBV-A3ab5ab was performed. Data points 
represent means with standard deviations of five eggs per condition, with all 
samples run and analyzed in triplicate using a ten-fold dilution series of IBV H52 BI 
as reference to determine virus quantity as EID 50 /ml equivalents [15]; (B) 
Embryonic death is indicated as a percentage of all remaining animals at each 
time point. 


3.5. rIBV-A3ab, rIBV-A3ab and rIBV-A3ab5ab are attenuated in vivo 
and provide protection against a homologous challenge 

To determine the in vivo phenotype of the generated recombi¬ 
nant viruses, one day-old SPF chickens were inoculated with via 
eye-drop with 10 3 EID 50 of IBV H52 BI, rIBV-wt, riBV-A3ab, rlBV- 
A5ab, or rIBV-A3ab5ab. No apparent clinical signs were observed 
in any of the birds during the course of the experiment, which 
was according to our experience upon inoculation under controlled 
conditions in isolators. As a readout for attenuation, the extent of 
ciliostasis in 10 tracheal sections per bird was therefore deter¬ 
mined at 7 days post inoculation. For rIBV-wt and IBV H52 Bi mean 



B 



Fig. 5. Ciliostasis after vaccination with IBV H52 BI and rIBVs, followed by challenge 
with IBV M41 in chickens. Plotted are tracheal ciliostasis scores per individual 
animal and means per experimental group. Maximal ciliostasis score per animal is 
40, which indicates complete ciliostasis in all 10 transversal tracheal sections 
examined. (A) Ciliostasis in one-day-old chickens determined seven days after 
vaccination with IBV H52 BI, rIBV-wt, rIBV-A3ab, rIBV-A5ab, rIBV-A3ab5ab, or 
non-vaccinated controls; (B) Ciliostasis in vaccinated and non-vaccinated animals 
determined seven days after challenge with IBV M41. A non-vaccinated non- 
challenged control group served as control. 
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ciliostasis scores were 26 and 19, respectively, while non- 
vaccinated control animals did not show any ciliostasis (Fig. 5A). 
Supplementary file 2 shows the ciliostasis induced by IBV H52 
wt and rIBV wt as observed in a comparable experiment (with 
scores of 34 and 40, respectively), demonstrating that both viruses 
are equally virulent for one-day-old chickens. All rIBV deletion 
mutants showed reduced ciliostasis compared to the wild type 
viruses, with rIBV-A3ab showing the most attenuated phenotype 
(Fig. 5A). A one-way ANOVA with HSD post hoc analysis showed 
that ciliostasis in IBV H52 Bl- and rlBV-wt vaccinated animals 
was significantly different from that in the non-vaccinated animals 
(P<.05), while ciliostasis scores from rIBV-A3ab, rlBV-A5ab, or 
rIBV-A3ab5ab vaccinated animals did not differ from those in 
non-vaccinated controls. These results indicate that the proteins 
encoded by the accessory genes 3 and 5 contribute to the 
pathogenicity of virulent IBV in young SPF chickens. 

At 21 days after vaccination, the remaining animals in each 
group were challenged via eye-drop with 10 3 E1D 50 IBV M41. Seven 
days post challenge, all remaining animals were euthanized, and 
their tracheal ciliary activity was evaluated to determine whether 
the chickens had been protected against the challenge with viru¬ 
lent IBV of the same serotype. The non-vaccinated animals showed 
a mean ciliostasis score of 39, indicative of a successful challenge. 
All vaccinated animals showed significantly reduced ciliostasis (p 
< .05 for all vaccinated groups compared to the group of non- 
vaccinated animals), with scores of 15 or lower for the recombi¬ 
nant IBV deletion mutants compared to 17 for IBV H52 BI. These 
data show that vaccination with recombinant IBV induced suffi¬ 
cient immune responses in the chickens to protect against a chal¬ 
lenge with virulent IBV of the same serotype. 


4. Discussion 

Targeted RNA recombination was performed to delete the 
accessory genes 3 and/or 5 from the IBV H52 Bl to thereby elimi¬ 
nate expression of the respective proteins. The intended recombi¬ 
nant IBVs were rescued; they propagated in ovo to titers similar 
to that of the parental virus. In chickens they exhibited reduced 
virulence and induced protective immune responses. The results 
demonstrate that our approach, which enables the genetic manip¬ 
ulation of a virus that typically cannot be propagated in vitro, pro¬ 
vides a new means to generate recombinant live attenuated 
vaccines against infectious bronchitis. 

Our RGS allows for selection of recombinant viruses based on 
their ability to propagate in embryonated eggs. After the first 
end-point dilution cycle, selection of the intended IBV recombi¬ 
nants was successfully accomplished based on IHC, RT-qPCR, and 
sequence analysis of the mutated regions. At that stage, we also 
observed for rIBV-A5ab and rIBV-A3ab5ab multiple amplicons 
for the mutated regions as well as spontaneous deletions in the 
3'-UTR (data not shown). Mutations and deletions in the hypervari¬ 
able region of the 3'-UTR, preceding a highly conserved region of 
293 nt at the 3' end of the 3'-UTR, have been observed in IBV field 
strains and laboratory strains earlier [30,31], but these viruses 
were not used for further studies here. 

Deleting the accessory genes from the IBV H52 genome did not 
affect the viability and replication capacity of the mutant viruses in 
ECEs relative to both wild type and recombinant wild type IBV 
H52. Apparently the proteins encoded by the accessory genes 3 
and 5 do not contribute critically to the replication of IBV in ovo, 
which is in accordance with previous observations in cell culture 
cells and in ECEs using recombinant Beaudette virus [9,16,17]. Pro¬ 
teins 3ab and 5ab have recently been associated with the chicken 
IFN response to IBV infection in cell culture [19-21], However, 
the ability to induce IFN as well as the sensitivity to the action of 


IFN is only fully functional during the third and last week of 
embryonic development [32,33] and develops in ECEs with gesta¬ 
tional age [33], This might explain the limited differences in 
growth kinetics of the rIBV deletion mutants in eight-day-old 
embryonated SPF chicken eggs. 

The proteins encoded by the accessory genes of virulent IBV 
contribute to the pathogenicity of the virus in chickens. Recombi¬ 
nant IBV viruses lacking the accessory genes have been generated 
before [9,16,17], but their contribution to virulence in chickens 
could not be addressed since the Beaudette strain used has no 
pathogenic phenotype in vivo 11,14]. The observed attenuation 
in chickens described here is in accordance with the observed phe¬ 
notypic characteristics of other coronaviruses in which these 
group-specific genes were inactivated. Deletion of the group- 
specific genes for mouse hepatitis virus MHV [23], feline infectious 
peritonitis virus FIPV [22], and transmissible gastroenteritis virus 
TGEV [24] resulted in reduction of viral loads, clinical symptoms 
and mortality after infection in their respective host (mouse, cat 
and pig, respectively). 

Coronaviruses that are attenuated by deletion of the accessory 
genes have been successfully used in vaccination trials. Thus, vac¬ 
cination with recombinant FIPV lacking the group-specific genes 3 
or, to a lesser extent, virus lacking genes 7, could protect cats 
against a lethal challenge with the parental virus [22], However, 
no protection was observed when both genes 3 and 7 had been 
deleted. In our study, rIBV-A3ab, rIBV-A5ab, and rIBV-A3ab5ab, 
each had the ability to protect against challenge with a virus of 
the same IBV serotype (Fig. 5B). In particular rlBV-A3ab showed 
the strongest attenuation based on the observed ciliostasis scores 
(Fig. 5A). A proper balance between attenuation and induction of 
protective immunity is critical when making a live attenuated 
virus vaccine and it needs to be tested in the appropriate animal 
model for each virus. In particular, the extent of attenuation caused 
by deletion of genes 3ab and 5ab might depend on the 
(pathogenicity of the) particular IBV strain being used. Finally, fur¬ 
ther studies will be needed to elucidate precisely the functions of 
the IBV 3ab and 5ab gene products in vivo, and how deletion of 
the genes results in the modulation of the host response to infec¬ 
tion and, eventually, to the reduced ciliostasis in chickens. 
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